APS/123-QED 



The Effect of the Pauli Exclusion Principle in the Many-Electron Wigner Function 

Emiliano Canccllieri, Paolo Bordone and Carlo Jacoboni 
Dipartimento di Fisica Universitd di Modena e Reggio Emilia and CNR-INFM 
S3 National Research Center, via Campi 213 /a, 1^1100 Modena (Italy). 



(Dated: February 6, 2008) 

An analysis of the Wigner function for identical particles is presented. Four situations have been 
considered, i) A scattering process between two indistinguishable electrons described by a minimum 
uncertainty wave packets showing the exchange and correlation hole in Wigner phase space, ii) An 
equilibrium ensemble of N electrons in a one-dimensional box and in a one-dimensional harmonic 
potential showing that the reduced single particle Wigner function as a function of the energy defined 
in the Wigner phase-space tends to a Fermi distribution, iii) The reduced one-particle transport- 
equation for the Wigner function in the case of interacting electrons showing the need for the two- 
particle reduced Wigner function within the BBGKY hierarchy scheme, iv) The electron-phonon 
interaction in the two-particle case showing co-participation of two electrons in the interaction with 
the phonon bath. 

PACS numbers: 05.30.Fk; 63.20.-e; 72.10.-d 



I. INTRODUCTION 

Highly sophisticated technologies produce physical sys- 
tems, and in particular semiconductor devices, of very 
small dimensions, comparable with electron wavelength 
or with electron coherence lengths. Under such condi- 
tions, semi-classical dynamics is not justified in principle 
and interference effects due to the linear superpositions 
of quantum states have to be considered. Among the 
possible different approaches, the Wigner- function (WF) 
has proved to be very useful for studying quantum elec- 
tron transport [H, 0, H, 01; owing to its strong analogy 
with the semiclassical picture, since it explicitly refers 
to variables defined in an (r, p) Wigner phase space, to- 
gether with a rigorous description of electron dynamics 
in quantum terms. 

In this work we present an analysis of the WF for iden- 
tical particles. Even thought the WF was defined from 
its very beginning for the study of many-particle physics, 
in electron transport theory it has been used mainly in its 
one-particle version. The importance of the many-body 
problem derives from the fact that any real physical sys- 
tem one can think of is composed of a set of interacting 
bodies. Moreover, since we are dealing with quantum me- 
chanical systems the symmetry properties that describe 
the behavior of identical particles play an essential role. 
The present paper will be focused mainly on the last sub- 
ject. 

In particular, four situations will be analyzed: i) A 
scattering process between two indistinguishable elec- 
trons described by minimum uncertainty wave packets, 
showing the exchange and correlation hole in Wigner 
phase space, ii) An equilibrium ensemble of N electrons 
in a box and in a harmonic potential, showing that the 
sum of the values of the WF that correspond to points 
in the Wigner phase-space with energy in a given inter- 
val, tends to a Fermi distribution, iii) The transport 



equation for interacting electrons, showing the BBGKY 
hierarchy when the integral, over the degrees of freedom 
of all the particles but one, are performed @,||. iv) The 
electron-phonon interaction in the case of two particles, 
where new Keldysh diagrams Q appear with respect to 
the one-electron case @ . 



II. WIGNER FUNCTION FOR MANY 
IDENTICAL PARTICLES 

The WF was introduced by Wigner in 1932 to study 
quantum corrections to classical statistical mechanics [l|, 
B EU EH ■ Thus, even though it is now used mainly in 
single particle problems, from the very beginning this 
function was defined for N particles as: 



fw(ri,Pi,-,r N ,p Ni t) = J rfsi...dsAre r^ s 'P* 

, / Si Sat \ 

XIP (^1 + —, ...,T N + ~^,tj 

xV> (ri- — ,...,rjv- -y,*J- 

(1) 

In the case of identical particles, the wave function de- 
scribing the many-body system satisfy well known sym- 
metry relations. When the position coordinates of two 
particles are interchanged the wave function remains un- 
affected (bosons) or changes sign (fermions). Since the 
WF is bilinear in the wave function it remains the same 
if the positions and, accordingly, the Wigner momenta of 
two particles are exchanged. 

This symmetry property of the WF allows the defi- 
nition of a reduced M-particle WF in a system of N 
particles as [IH, EH : 
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fw ( r i>Pi>- 



, I'M, PM>t) — 



(N - M)lh 3 ( N - M ) 



dr M +i dp M+1 ...dr N dp N f w (r 1 ,p 1 , ...,r N ,p N ,t), (2) 
I 



where the superscript (N) indicates that the reduced M- 
particle WF is defined in a system with N particles. Note 
that in the case where M = 1 the above equation be- 
comes: 



fw ( r i>Pi»*) 



N 



dv2 dp2---dYM dp 



h 3(N-l) 

x/w(ri,Pi, ■■■,r N ,p N ,t). 



N 



(3) 



The factorials appearing in front of the integral in equa- 
tion ([2]) simplify to N in equation ([3]) since this is the 
number of equivalent ways one can reduce the TV-particle 
WF when the particles themself are supposed to be iden- 
tical. 



The WF for Many Single-Particle Wave 
Functions 



where the upper sign is for bosons and the lower for 
fermions. In the WF expression it is possible to identify 
two different types of terms. The first one is character- 
ized by the product of single-particle WFs. In each of 
these contributions, from the different wave functions, 
N WF are obtained that are evaluated in a particu- 
lar permutation of the variable indices as, for example: 
fw! (i"4 , P4 ) fw 2 (r i , Pi ) fw 3 (r 2 , P2 ) ■ ■ ■ fw N (r N -5 , Pns) ■ 



We consider the case of N particles in the system and 
we define the WF with a wave function that is a symmet- 
ric or anti-symmetric linear combination of single-particle 
wave functions V*( r )) (* = 1> — > ^0 as 



t/>(ri,...,rjv) = Vi(ri)V'2(r2)V'3(i"3) ■■•i/'Jv( r Af) 



±V>i(r 2 )V>2(ri)V'3(r3) -^n{pn) 
+f/>i(r 2 )?/'2(r3)V'3(ri) ...ip N (r N ) 



(4) 



The second type of contributions accounts for the ex- 
change effects and vanishes when the wave functions 
ip n ( r ) do not overlap. These terms are constituted by in- 
tegrals of the product of N factors tp n one for each of 
the TV wavefunctions tp n . In these terms at list two prod- 
ucts ipn{^i + s i/2)V ; n( r i — s i/2) are evaluated with i ^ j. 
It is the presence of such factors that makes impossible to 
obtain the many-particle WF in terms of single-particle 
WFs. The number of factors ^ n f/>*, where ip n and if>* 
correspond to different particles, appearing in a given in- 
tegral can range from 2 to N. As an example the WF in 
the case of N — 2 reads: 



JW(ri,Pi,r 2 ,P2,) 



•/Vi(ri,Pi)jV 2 (r 2 ,P2) + /wi(r2,p 2 )/w2(ri,pi) 
±JL J d Sl ds 2 e-^ SlPl+S2P ^ 

, l(ri + |)«( r2 -|), 2 ( r2 + |),.( ri -|) 



(5) 



here 4 terms appear, 2 for each kind of contribution. The 
two-particle system is treated in details in [14j |. 



B. Example of Two Colliding Electrons 

A one-dimensional situation where two fermions col- 
lide with each other has been simulated. The Schrodinger 
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FIG. 1: One dimensional reduced single-particle WF of two 
interacting electrons at different times. The figure clearly 
shows the exchange hole due to the Pauli exclusion principle. 

equation was solved with initial conditions given by two 
minimum-uncertainty wave packets interacting through 
the Coulomb potential and the WF was evaluated at dif- 
ferent time steps (see Fig. [I]). 

In this figure we plot the one-particle reduced WF of 
the system for the case of two Gaussian wave packets with 
opposite central wave vectors. Since we are dealing with a 
one dimensional system, the two particles are expected to 



decelerate, scatter, and then move away from each other. 
At t = we suppose the two particles to be described 
by an antisymmetric wave function. In Fig. l.B the 
system is shown 2 ps after the Coulomb interaction is 
switched on. At the beginning of the scattering process 
the exchange hole due to the Pauli's exclusion principle 
appears. In part \.C the two particles are shown when 
their mutual distance has reached the minimum value. In 
this case the exchange hole is maximally evident. When 
the two particles are moving far enough from each other 
the exchange hole tends to disappear (part l.D). 



III. EQUILIBRIUM WF FOR NON 
INTERACTING PARTICLES IN CONFINED 
POTENTIALS 

In this section a system of N fermions in a confined 
potential has been studied. In the non-interacting case 
the one-particle reduced WF is studied at thermal equi- 
librium at a temperature of T = 2 K. In order to simplfy 
the mathematical treatment we shall introduce the sec- 
ond quantization notation. The ./V-particles wave func- 
tion can thus be written as: 



V(r 1; ...,r N ) = (r l7 ...,r N \iP) = (0|$(n)...¥(r w )M, 

(6) 

and the WF as: 



N 



h 3(N- 

(o|# (ri + ^) (r N + 



SJV 
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SjV 



ri 



si 



(7) 



where, here and in the following, ^> and W are the cre- 
ation and annihilation field operators. Since we are inter- 
ested in the thermal equilibrium distribution of a fixed 
number of particles, the density matrix in the above equa- 
tion is: 

where Z is the partition function, H the Hamiltonian, 
ks the Boltzman constant, and T the temperature of the 
system. Since the particles in the system are supposed 
to be non interacting, and the system is supposed to be 



confined, the Hamiltonian in its second quantization form 
can easely be written in terms of the particles creation 
& n and annihilation operators c„ as: 



H = ^2 e n clc n , (9) 

n=l 

where e is the energy of the n-th discrete level. Writing 
the field operators in terms of the creation and annihi- 
lation operators, the mean value element appearing in 
equation ([7]) becames: 
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0) x 



^ ( ri + y) + ^ ( rN 



S N 



41, („ - S {) 



0)e 



kg 1 •• n 1 n N > 



n\ n'„ nv 



^ (n + |) ...^ (r* + f ) ^ (r w - f ) (n - I) 



(10) 



where indicates the n-th eigenstate of the confining 
potential. In order to better understand how to treat 
the above equation let us focus on the two-particle case. 
Using the fermionic or the bosonic commutation rules for 
the contribution containing the creation and annihilation 
operators, the following identity is achieved: 



By substituting equation (fTT|) into equation (flQ| for 
N = 2, and using it in equation ([7]) the following expres- 
sion is obtained: 



(^™i n i'^ n 2 n 2 ^ ^™i n 2"^™2™'l' 



(11) 



/^(n,Pi) = j^J dr 2 dp 2 J ds 1 ds 2 e~ ± n^ +s ^Y.Y. e ~ 1 ^ {en ' 1+€ "' 2) 

Th\ Tl' 2 

± A, (r, + |) *„; (r, + «, (r, - |) C; (« " |) } , 
then, performing the integrals, we finally get: 



/#(ri,Pi) 



|EE e ^ (e "' i+x) /^ (ri ' Pi) 



4ee< 



; R i V 71.. 7l ' 



| ds ie -i""i. n . (r, + |) «, (r, - 5. 



e -b "i "2' ± e 



. 1 ^ 2e / 



JW „/ (ri,Pi), 



(13) 



n I n~ 



where /v^„ ( r iiPi) indicates the WF of the n\ h eigen- 
state of the confined potential. 

This expression can be written in a more compact form 
by identifying the bosonic and the fermionic cases: 

«1 "2 

(14) 



/$_ io „,KpO = !E E e-'^/^Kpi). 

(15) 

The generalization to the iV-particle system is straight- 
forward and gives: 
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, + £„,+.. .+e 



(17) 



Let us note that in the above expression for the 
fermionic case the following sum appears: 



'T T T T | — T — |T 



• « 4 particles 

— 6 particles 
" 8 particles 
»— x 10 particles 

— Fermi function 




20 

Energy K 



40 



, +«»„+.. •+« 



(18) 



In the limit of large N and of an infinite number of al- 
lowed states with continuous energies the above term 
gives the Fermi function evaluated at an energy value 
e ni PI. 



FIG. 2: Average values of the points corresponding to the 
same energy inteval of the single-particle reduced WF of a 
system of N particles at thermal equilibrium at a temperature 
of 2 K in an ID infinite square well potential. Since the energy 
depends only upon the momenta of the particles, the above 
average corresponds to the integral over the position variable 
(x). The width of the well has been fixed to 150 nm. In 
the upper part of the figure the black triangles indicate the 
energies corresponding to the eigenstates of the well. 



A. Infinite Square Potential Well 



An inifite square well potential in one dimension has 
been investigated at a temperature of T = 2 K. The 
single-particle WF has been evaluated for N = 4,6,8 
and 10 by means of equation (|17p. Then the average 
values of the points of the WF corresponding to energy 
interval e, e + Se have been plotted as a function of e, in 
Fig. [21 Since the energy depends only upon the Wigner 
momentum of the particle (p 2 /2m), the above average 
corresponds to the integral of the single-particle reduced 
WF with respect to the position variable (x). In our 
simulations the width of the well has been kept contant 
to a value of 150 nm. 

A comparison between our curves and the Fermi func- 
tions is obtained by evaluating the chemical potentials fj, 
for N — 4, 6, 8 and 10 from a numerical solution of the 
equation: 

oo 

^ = E^Z5 ■ ( 19 ) 

n =l e k B T + 1 

The curves in Fig. [5] show a very good agrement between 
the Fermi function and the average of the WF for any 
number of particles. 

It is worth noting that, as expected, the agreement 
between the averages of WFs and the Fermi distributions 
is higher as the number of particles increases. However 
in Fig. [5] even in the case of 10 particles the value of the 
WF's average corresponding to the point with e = do 
not reach the maximal value of 1 . For this reason we have 



0,98- 




Energy K 



FIG. 3: Average values of the points corresponding to the 
same energy of the single-particle reduced WF of a system of 
N particles at thermal equilibrium at a temperature of 2 K in 
a ID infinite square well potential. The width of the well has 
been fixed to 150 nm and the number of particles increased 
from 10 to 85. 



plotted in Fig. [3] the system in more ditail in the energy 
range from to 10 K for a higher number of particles. As 
before we simulate a well 150 nm wide with an electron 
gas at a temperature of 2 K, in this case, however, the 
number of simulated particles is increased to 85. Fig. [3] 
shows that when the number of particles increases the 
value corresponding to e = approaches 1 . 
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Energy K 

FIG. 4: Average values of the points corresponding to the 
same energy of the reduced single-particle WF in a harmonic 
potential with spring constant k = 1.54 x 10 s [kg/sec 2 ]. The 
system in the case of N = 4, 6, 8, 10 particles is studied at 
thermal equilibrium at a temperature of 2 K. The curves 
clearly tend to the Fermi-Dirac distribution. In the upper 
part of the figure the black triangles indicate the energies 
corresponding to the eigenstates of the harmonic potential. 



°0 20 40 60 80 100 120 " 140 " 160 180 
Energy K 

FIG. 5: Average values of the points corresponding to the 
same energy of the reduced single-particle WF in a ID infi- 
nite square potential well. In a system with 10 particles at a 
temperature of 2 K, the width of the well has been reduced 
from 150 nm to 70 nm, that corresponds to an energy gap 
decrease, from the ninth to the tenth energy level, from 3.8 to 
17.1 K. When the energy spacing between the levels is bigger 
than the thermal energy, the particles distribution deviates 
form a Fermi function. 



B. Harmonic Potential 

As a second example, we have studied a one- 
dimensional harmonic potential. Equation p7[) has been 
evaluated for different numbers of particles (4, 6, 8 and 
10) at a temperture of T = 2 K. The average over the 
points of the WF belonging to the same energy interval 
(p 2 /2m+^kx 2 , where m is the mass of any particle in the 
system) are plotted in Fig |U The Fermi function with 
the chemical potential given by equation (|19p is clearly 
approached by the corresponding average of the WF. 



to 2 K. The simulations show that such oscillations get 
more and more evident as the ratio between the spacing 
of the energy levels and fcsT becomes greater. Under 
these conditions the picture of a continuous spectrum of 
energies breaks down. 

IV. TRANSPORT EQUATION 

The dynamical equation for the single-particle WF is 
derived by differentiating the definition of the WF itself: 



C. Effect of energy separation between levels 

It is possible to study how the particle distributions 
change when the dimension of the well or the strength 
of the harmonic potential are varied in the case of both 
the infinite square well potential and the harmonic po- 
tential. When the width of the well decreases or the 
strength of the harmonic potential increases, the spacing 
between the allowed energy levels increases and an oscil- 
lating behaviour shows up in the curves (see Fig.sJSl [HI- 
Our calculations have been performed in the case of a 
system with 10 particles at a temperature of 2 K. When 
the width of the well is decreased from 150 nm to 70 
nm, the energy gap between the ninth and the tenth en- 
ergy level increases from 3.8 to 17.1 K. In the case of the 
harmonic potential the strength of the force constant is 
varied from 1.54 x 10 s to 5.69 x 10 s kg/sec 2 leading to an 
increase of the distance between the energy levels from 1 



. ( 2 °) 

By means of the Schrodingcr equation it is possible to 
evaluate the time derivative of the product of the two 
wave functions and to obtain the dynamical equation for 
the WF Q: 

d p 

—fw{r,p,t) = VjV(r,p,i) 

at m 

+ 7^ J dp'V w (r,p-p')f w (r,p',t), 

(21) 

where Vw is the interaction kernel for an external poten- 
tial V(y). Note that the interaction term, given by: 



W(r,p)=i / </s< 



(22) 
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FIG. 6: Average values of the points corresponding to the 
same energy interval of the reduced single-particle WF in a 
harmonic potential. In the case of 10 particles the bound 
constant has been varied from k — 1.54 x 10 8 [kg/sec 2 ] to 
k = 5.69 x 10 8 [kg /sec 2 ] corresponding to an increase of the 
spacing between the energy levels from 1 to 2 K. When the gap 
between the energy levels increases and becomes bigger than 
2 K, the gas temperature, the particles distribution shows an 
oscillating behaviuor superimposed to the Fermi-like shape. 



depends on the values of V at points different from r. 
However, while the non-locality of Vw extends to infin- 
ity, its effect on the electron dynamics has to be consid- 
ered only up to regions where the electron correlation is 
different from zero. 



A. Electron-Electron Scattering 



Let us study the transport equation for electron- 
electron scattering. In the case where no phonons nor 
external forces are present, the potential V(ri,r2...rjy) 
is given by the Coulomb interaction, the transport equa- 
tion reads: 



7j7/w(ri,pi, ...,r N ,p N ,t) 



Y] — V r Jw(ri,Pi, -,ri,pi, ...,r N ,p N ,t) 
* — ' m 

I 

ffYH2j d Pi d Pj s (. A Pi+ A Pj)Vw(\Ti-Tj\,&Pi-&Pi) 



i 3 



x/w(ri,pi, p-, ...,r,-,p' , ...r N ,p N ,t), 



(23) 



where Vw is the potential kernel of the Wigner equation 
and Ap = p p'. As done before, in order to get a 
better understanding of the above equation, the kernel 
that describes the electron-electron interaction is studied 
for N = 2: 



/ (r 1 ,r 2 ,p 1 ,p 2 ) = — J d Sl ds 2 e i^" lT '«' 

xV {v 1 + -,v 2 + -)-V{r 1 --,v 2 --) . 

(24) 



Since the Coulomb interaction depends only upon the 
distance between the two particles it is usefull to re-write 
the above equation using the new variables x = ri — r 2 , 
s = si - s 2 , and s' = (pis 2 4- p 2 si)/(pi + p 2 ): 
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W(ri,r 2 , Pl ,p 2 ) = 1 J dBds / e-*( pi+Pa )"'e-t<P 1 - pa >[v(x + |) -v(x 



i 

^_ 

ih 

h 3 5(pi + p 2 )VV(n - r 2 ,pi - p 2 ). 



5(Pi+P 2 ) / dse-*(P 1 -P a >[7(x+|) -^(x-f) 



(25) 



Thus the factor <$(A P j + A Pj ) appearing in equation 
(|23|) represents the constrain for the total momentum 
consevation while the difference (Ap, — A P j) indicates 
that the interaction depends only upon the momentum 
transfer between particle i and j. 

When the single-particle reduced WF in the case of N 
particles is evaluated, equation (f2"5)) reads: 

^(r.p.t) = -£v r /™(r lPl t) 

+J_ J dgdp e J dp'V w (\r - g\,2Ap) 

X/w( r 'P')^Pfi'*)) ( 26 ) 

where r and p are the position and the Wigner momen- 
tum of the considered particle and, g and p g indicate 
the position coordinates of one of the remaining N — 1 



particles. It should be noticed that all the particles are 
interacting with each other, but, due to their indistin- 
guishability, all the contributions are identical and sum 
up to balance the factorials appearing in equation ((2). 
The above expression shows that the transport equation 
for the reduced single-particle WF depends on the re- 
duced two-particle WF. When the transport equation for 
the reduced two-particle WF is evaluated, the electron- 
electron interaction term depends upon the three-particle 
reduced WF and so on for the transport equation for 
the other reduced WFs. It is the Wigner picture of the 
BBGKY hierarchy. 

When the WF is written in terms of antisimmctric 
single-particle wave functions, as we have seen in equa- 
tion |5]). two types of contributions can be identified. It is 
then possible to study how the transport equation reads 
when only the contributions due to non overlapping wave 
functions are considered: 



d_AN) 
dt Jw 



f^(r,p,t) = ^£.C(r,P,*) 



= -£v r /#°(r,p,*) + ^£ / dp'dgV w (\r-g\,2Ap) 



x/g?(r,p',t). 



(27) 



Besides a Liouvillian contribution, an interaction term B. Electron-Phonon Scattering for the 2-Electrons 

appears where each one-particle contribution interacts WF 
with all the others as in the Hartree approximation. In 
the case of overlapping wave functions also the other kind 

of contributions (as studied in equation ((SJ)) must be con- The e-ph interaction term for the dynamical equation 

sidered, and the exchange term is restored. for two electrons is made out of eight terms as follows: 
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di fw 



ep q' 

7w fr x ,Pi + ^-,r 2 ,P2,{-,n q ' - 1, ...},{n' q },t 
+e -i(q'n-^(t-to))^ ( + 1/w ^ ri)Pl _ ^, r2 , pa, {„,},{„., + 1, ...},* 



i(q'n-« /(i-t )) /Z7 



"g' /w fri,pi + ^,r 2 ,P2,{n q },{...,n' q , -l,...},tj + o.p.j, 



(28) 



r 



where ri,pi,r2,p2 are the Wigner phase space coordi- 
nates of the two particles. In the above equation o.p. 
stands for other particle and indicates the four terms 
where r 2 replaces ri in the exponential factors and p 2 un- 
dergoes a variation of ftq'/2 while pi remains unchanged. 

The eight terms appearing on the r.h.s. of the above 
equation have simple physical interpretations: the e-ph 
interaction occurs as emission or absorpion of a quantum 
of any mode q and this may appear in the state on the 
left or on the right of the bilinear expression that defines 
the WF. Each elementary interaction or vertex changes 
only one of the two sets of variables of the WF; more 
precisely, one of the occupation numbers n q is changed 
by unity and one of the electron momenta is changed by 
half of the phonon momentum. 

In analogy with the Chambers formulation [16[ of the 
classical kinetic equation it is possible to introduce new 
variables (r*,p*,t*) that allow us to obtain an integral 
form of the dynamical equation for the WF. This inte- 
gral equation is in a closed form and can be solved by 
iteratively substituting it into itself, leading to what is 
known as its Neumann expansion. 

Equation (|2"5)) gives 8 terms for the contribution of 
the first order of the Neumann expansion, 64 terms for 
the contribution of the second order and so on for the 
higher order terms. In order to obtain meaningful physi- 
cal quantities, however, the trace over the phonon modes 
must be performed, leading to a vanishing contribution 
for each term corresponding to an odd order in the Neu- 
mann expansion [8[ . Only terms with an even number of 
verteces give dyagonal (in the phonon modes) contribu- 
tions different from zero. 

As stated before, the second order in the Neumann 
expansion gives 64 terms. Among these, 32 yield con- 
tributions dyagonal in the phonon modes and survive to 
the trace operation, 16 terms refer to one particle and 16 
to the other. For each particle 8 terms are the complex 
conjugate of the other 8 and can be summed together 
leading to 8 contributions for each particle. Among these 
it is possible to recognize four standard interactions un- 



dergone by each particle: real emission; real absorption; 
virtual emission and virtual absorption. 

The main difference to the single particle case lies in 
the eight (four for each particle) remaining terms. In 
the two-particle case the phonon occupation number can 
be changed, in the first or second set of values in the 
arguments of the WF, not just by one electron loosing 
(gaining) a Wigner momentum equal to half the phonon 
momentum in each of the two verteces Q but also by the 
action of two electrons. One electron looses or gains half 
phonon momentum in the first vertex and another elec- 
tron looses or gains half phonon momentum in the second 
vertex. Since we are dealing with identical particles we 
don't know which electron interact with the phonon bath 
in the first or in the second vertex. These four terms are 
real or virtual emissions or absorptions where the inter- 
action with the phonon bath is shared between the two 
electrons. 

It should be recalled that the p variable of the WF 
is obtained as a linear combination of two electron mo- 
menta. A specific value p is obtained as p = ft (ki +k2)/2 
where ki and k2 range from — oo to +oo. For this reason 
the Wigner momentum undergoes a change correspond- 
ing to half of the phonon momentum at each interaction 
vertex. 

In Fig. by means of the Keldysh-diagram formal- 
ism, two of the co-participated graphs are shown: a 
real phonon emission due to electron-electron coopera- 
tion and a virtual emission respectively. In each case both 
the Keldysh diagram representing the transition and the 
corresponding Wigner path [8| undergone by the two elec- 
trons are shown. Since the Keldysh diagrams are in the 
density matrix (DM) representation, four time-lines ap- 
pear, two for each electron. The first and the third lines 
correspond to the first wave function of the DM while the 
other two correspond to the second wave function. In the 
upper box at time t = t\ one electron in the first wave 
function of the DM emits a phonon and at time t = t 2 
another electron in the second wave function emits the 
same phonon. 
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In the lower box of Fig. [7] at time t = ti an electron 
absorbs the phonon emitted at t = ti, corresponding to 
a virtual emission. 



V. CONCLUSION 

We have developed a model based on the WF formal- 
ism that allows to introduce the symmetry effect in a sys- 
tem where electrons interact with each other and with the 
phonon bath. We have shown how this formalism can be 
uscfull by applying it to different situations: the study of 
the electron-electron scattering, the study of the thermal 
distribution of N particles in confining potentials, and 
the study of the two-electron dynamics in the presence 
of electron-phonon scattering. We have also shown that 
with the WF it is possible to reproduce a Fermi like dis- 
tribution defining the energy in the Wigner phase-space. 



FIG. 7: Real emission of a phonon mode in mutual partecipa- 
tion by two electrons. One electron ei changes its Wigner mo- 
mentum by half of the phonon momentum at time t = t\ and 
the other electron e2 does the same at a later time t = ti (up- 
per bow). Virtual emission where one electron lose a Wigner 
momentum equal to half the momentum of the phonon while 
the other electron gains the same amount (lower bow). 
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